Thin helium film on a glass substrate 
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We investigate by Monte Carlo simulations the structure, energetics and superfluid properties 
of thin 4 He films (up to four layers) on a glass substrate, at low temperature. The first adsorbed 
layer is found to be solid and "inert", i.e., atoms are localized and do not participate to quantum 
exchanges. Additional layers are liquid, with no clear layer separation above the second one. It is 
found that a single 3 He impurity resides on the outmost layer, not significantly further away from 
the substrate than 4 He atoms on the same layer. 

PACS numbers: 67.25.dp, 67.60.-g, 02.70.Ss 



INTRODUCTION 

Adsorption of thin films of the most abundant isotope 
of He ( 4 He) on various substrates, is a subject of long- 
standing interest, motivated by the unusual, often in- 
triguing physical properties that these films display [![. 
A chief example is the rich phase diagram of Helium on 
graphite, which has been the subject of an intense investi- 
gation over the past three decades. Significant attention 
has also been devoted to adsorption on weakly attractive 
substrates, e.g., those of alkali metals, following origi- 
nal suggestions that, at low temperature, the weakness 
of the physisorption potential ought to result in unusual 
wetting [3, 0| and superfluid [3] behavior. 

Less is known quantitatively, both theoretically and 
experimentally, about adsorption over a glass substrate, 
despite its obvious relevance to experimental studies of 
the superfluid properties of liquid (and solid) helium em- 
bedded in such porous materials such as vycor or aerogel 
[H . Of particular interest are issues such as the interplay 
between the one or more solid, "inert" layers that are 
expected to form over the substrate, and the liquid-like 
behaviour that should progressively emerge as successive 
layers are adsorbed. Glass is the archetypal disordered 
substrate, and that renders its theoretical modelling par- 
ticularly challenging. 

Microscopic theoretical studies of adsorption of helium 
(both pure 0, 0| and in isotopic mixtures on glass, 
based on realistic interaction potentials, have typically 
made use of a vary simple model of the substrate, de- 
scribed as a flat surface on account of its lack of a definite 
crystal structure. As mentioned above, the substrate is 
in reality irregular, disordered, and this can be expected 
to promote the appearance of a stable, low-coverage in- 
sulating film, in which atoms are localized at randomly 
located wells of the adsorption potential; such a phase is 
obviously lost if a flat substrate model is assumed. As 
the coverage is increased, however, a first-order transi- 
tion to a regular (monolayer) crystal is expected, and as 
further layers are adsorbed the corrugation of the sub- 
strate ought to become less relevant, i.e., a flat substrate 
model should capture the essential physics of the prob- 
lem (this is what observed, e.g., in numerical simulations 
of helium films on regular crystalline substrates, such as 



graphite [§]). 

In this work we study the adsorption of 4 He on a 
glass substrate, modelled as in all previous studies by 
others, by means of numerical simulations based on the 
continuous-space Worm Algorithm (WA). The purpose 
of this investigation is that of determining quantitatively 
structural and superfluid properties of the adsorbed 4 He 
film at low temperature. Specifically, we present results 
for T=l K and T—0.5 K. We considered films up to four 
layers thick, and also looked at a single 3 He impurity di- 
luted in the film, for different values of coverage. There 
is considerable uncertainty as to the actual values of the 
parameters that are used in the potential describing the 
interaction of a helium atom with the substrate, and we 
attempted to pick values that overlap with most previous 
studies. 

Our main findings are the following: 

1. The first helium adlayer is solid, whereas the suc- 
cessive ones are liquid. In the solid "inert" layer, 
atoms are highly localized, the layer is fairly well 
separated from the second, and quantum exchanges 
(both intra- and inter-layer) are virtually absent. 
On the other hand, successive layers are liquid- 
like, with no sharp separation between layers above 
the second. Atomic exchanges take place, mostly 
within the layer for the second layer, as well as 
inter-layer above the second one. Layers above the 
first are all superfluid in the low temperature limit. 
Thus, one can conceptually distinguish the first ad- 
layer from the rest of the system. Previous studies 
[8(] have often assumed that two inert layers would 
form, but the results presented here show that the 
model potential utilized is not sufficiently strong to 
stabilize more than one such layer. 

2. A single 3 He impurity diluted in the film, for differ- 
ent values of 4 He coverage, always positions itself 
on the top layer. Its wave function is not signifi- 
cantly more spread out than that of a 4 He atom on 
the same layer, i.e., the 3 He atom does not spend 
a significantly fraction of time at greater distances 
from the substrate. In other words, 3 He atoms do 
not "float" on top of the 4 He film, as one might 
have expected. 
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The remainder of this manuscript is organized as fol- 
lows: in the next section, the model Hamiltonian is 
briefly introduced, and some details of the computational 
method utilized are presented (for a thorough illustra- 
tion of the methodology, the reader is referred to Refs. 
loL 11); we then discuss our results and outline our con- 
clusions in the two following sections. 



MODEL 

Our system of interest is modelled as an ensemble of 
N 4 He atoms, regarded as point particles, moving in the 
presence of an infinite, smooth planar substrate (posi- 
tioned at z=0). The system is enclosed in a vessel shaped 
as a parallelepiped, with periodic boundary conditions in 
all directions [12J. Let A be the area of the substrate; 
correspondingly, the nominal 4 He coverage is 9 = N/A. 
The quantum-mechanical many-body Hamiltonian is the 
following: 



N 



N 



(i) 



Here, A=6.0596 KA 2 , V is the potential describing the in- 
teraction between two helium atoms, only depending on 
their relative distance, whereas U is the potential describ- 
ing the interaction of a helium atom with the substrate, 
also depending only on the distance of the atom from 
the substrate. We use the accepted Aziz potential 13 1 
to describe the interaction of two Helium atoms, which 
has been used in almost all previous studies, and has also 
been shown 11| to provide an accurate description of the 
energetics and structural properties of liquid 4 He in the 
superfluid phase. 

More delicate is the issue of which potential to use 
to describe the interaction of a helium atom with the 
smooth substrate (i.e., the U term in (JTJ)). The simplest 
potential is the so-called "3-9" : 



D fa 



,9 



(2) 



which is a functional form obtained by integrating the 
Lennard- Jones potential over a semi-infinite, continuous 
slab. The parameters a and D are normally adjusted 
to fit the results of some ab initio electronic structure 
calculations for the specific adatom-substrate system of 
interest. D is the characteristic depth of the attractive 
well of the potential, whereas a is essentially the location 
of the minimum of such a well. 

There appears to be some ambiguity as to the values 
of the coefficients that would make this potential appro- 
priate for the system of interest here. Choices made in 
some of the previous studies are summarized in Table Q] 
Although these studies were aimed at describing the en- 
vironment experienced by a helium atom in the confines 
of various porous materials, such as vycor, the "build- 
ing block" for all of them is the interaction potential felt 



TABLE I: Values of the coefficients a and D that appear in the 
"3-9" potential describing the interaction between a helium 
atom and a glass substrate, utilized in previous studies. 



Ref. 


a (A) D (K) 


6 


2.18 128 


7 


3.6 87 


8 


2.05 100 



100 




FIG. 1: Comparison of the helium-glass adsorption potential 
as given by expression [2] with the parameters provided in Ref. 
H| (solid line) and0 (dashed line). 



by a helium atom close to an infinite, flat glass (SiC^) 
substrate. As is seen in the table, there seems to be a rea- 
sonable agreement over the value of D, which is ~ 100 K 
in all studies. On the other hand, while a has been taken 
to be slightly above 2 A in Refs. @, [H, in Ref. it was 
set to a much larger value, namely 3.6 A. As shown in 
Fig. [TJ the ensuing potentials are quantitatively very dif- 
ferent, the parameters utilized in Ref. 7; giving rise to a 
potential with a significantly more attractive long-range 
tail (in turn likely to give rise to quantitatively different 
adsorption characteristics) . 

The claim is made in Ref. that a=3.6 A applies to 
the interaction of helium atoms with a substrate of MgO 
(regarded as sufficiently similar to Si02); but the most 
reliable estimates reported in the literature (see, e.g., Ref. 
\l4 ) place a for such a system in the 2.2-2.5 A range. In 
this work, we utilize the same values of a and D taken in 
Ref. [H, which are close to those of Ref. [H 



METHODOLOGY 

We performed Quantum Monte Carlo (QMC) simu- 
lations of the system described by Eq. (TT]) using the 
continuous-space Worm Algorithm. This methodology 
has emerged over the past few years as a far superior 
alternative to conventional Path Integral Monte Carlo 
(PIMC), to compute accurate thermodynamic properties 
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FIG. 2: Energy per particle in the low-temperature limit for 
an adsorbed 4 He film on a glass substrate, computed by QMC 
simulations (circles), as a function of coverage. Statistical er- 
rors are smaller than the symbol size. Dashed line through 
circles is a polynomial fit to the data. Dotted line is the chem- 
ical potential, obtained from the fit. The label S+L indicates 
a region of liquid puddles on top of the first solid layer. 
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FIG. 3: Density profiles n(z) in A -3 , computed by QMC, 
in the direction perpendicular to the substrate, at various 
coverages. The temperature is T=l K. 



Energetics 



of Bosc systems at finite temperature. Specifically, the 
WA allows one to simulate systems comprising up to two 
orders of magnitude more particles than PIMC, allowing 
for a reliable extrapolation of values of physical observ- 
ables to the thermodynamic limit. 

The specific implementation utilized in this project is 
canonical, i.e., we keep the number N of particles fixed 
[lit ■ Other technical aspects of the calculations are com- 
mon to any other QMC simulation scheme. 



RESULTS 
Computational details 

Results were obtained by simulating systems compris- 
ing a number N of particles ranging from 36 to 108, span- 
ning a range of 4 He coverage 9 from 0.050 to 0.300 A~ 2 . 
Estimates presented here are obtained at two tempera- 
tures, namely T=l K and T=0.5 K; based on the wealth 
of accumulated theoretical and experimental data for liq- 
uid 4 He, T=l K is low enough that structural and ener- 
getic properties are essentially those of the ground state. 
The usual fourth-order high-temperature propagator uti- 
lized in all previous studies based on the WA was adopted 
here; convergence of the estimates was observed for a 
time step r w 3.125 x 10~ 3 K _1 , which corresponds to 
P=320 imaginary time "slices" at T=l K. Because the 
computational cost was negligible, all estimates reported 
here were obtained using twice as small a time step, in 
order to be on the safe side. 



Fig. [J] shows the energy per particle e(9) as a function 
of coverage, computed by QMC simulations at T=l K. 
Assuming these to be essentially ground state results, 
based on a piecewise polynomial fit to the estimates for 
e(9) and the relation 



K0) = e(9) + 9- (3) 



we obtained the T=0 chemical potential, also shown in 

Fig. rj 

The monolayer equilibrium coverage 9 e , for which 
fj,(9 e ) — e(9 e ), is found to be close to 0.052 A -2 , to be 
compared to 0.043 A~ 2 for the purely two-dimensional 
system This corresponds to a superfluid monolayer, 
but it is important to restate that the model of substrate 
considered here lacks the corrugation and disorder of the 
surface, which will stabilize a disordered, "glassy" insu- 
lating monolayer in the real system. On increasing cov- 
erage, the system forms a two-dimensional (triangular) 
crystal at 9\c ~ 0.08 A~ 2 . The crystal can be com- 
pressed up to 6>ip=0.100 A~ 2 , at which second layer pro- 
motion is observed. The system forms a stable liquid 
second layer at 6> 2 l=0.153 A~ 2 , whereas in the coverage 
region 9\p < 9 < 9il puddles of liquid form on top of 
the solid first layer. 

Above 9il , the chemical potential is a smooth, monotoni- 
cally increasing function of coverage, approaching asymp- 
totically the value ~ -7.2 K for bulk liquid 4 He. This 
signals continuous growth of the adsorbed film, without 
well-defined layer formation. This is confirmed by direct 
observation of the structure of the film, as well as by its 
superfluid properties. 
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FIG. 4: Integrated coverage 00) (in A" 2 ), defined in Eq. H 
as a function of the distance z from the substrate, for a film 
of total coverage 0.299 A -2 . The temperature is T=l K. 




e (A 2 ) 

FIG. 5: He superfluid fraction ps of adsorbed films of various 
coverages 6 (in A -2 ), computed by QMC simulations. The 
temperature is T=0.5 K. 



Structure 

Fig. [3] shows the computed 4 He density profiles n{z), 
in the direction perpendicular to the substrate, for four 
different coverages. The first observation is that there 
exists a fairly sharp demarcation between the innermost 
(first) layer and the second, with essentially no overlap 
between the two, and no interlayer atomic exchange. On 
the other hand, second and third layers already overlap 
significantly, and essentially no separation exists between 
third and fourth layers. The fourth layer is broad and 
ill-defined, its average density very close to that of equi- 
librium liquid 4 He (0.022 A" 3 ). 

Additional evidence of the ill-defined character of lay- 
ers in this system is provided in Fig. |4j where the inte- 
grated coverage 

6{z) = J dx' dy' J dz'p(x',y',z') (4) 

which depends on the distance from the substrate, is 
shown for a film of total coverage 61=0.299 A" 2 . Steps 
can be seen in correspondence of the appearance of the 
second layer (at approximately 0.100 A -2 ) and of the 
third at roug hly 0.165 A~ 2 can be clearly seen. How- 
ever, the boundary between third and fourth layer is al- 
ready scarcely visible. This is not qualitatively dissimilar 
from what is observed on other substrates [12]. or The 
wealth of theoretical data accumulated over a wide va- 
riety of substrates, generally suggests that although the 
details of the first layer clearly reflect the properties of the 
underlying substrate (strength and corrugation), beyond 
the first layer what one observes seems to be essentially 
substrate-independent . 

Atoms in the first layer are arranged on a regular tri- 
angular crystal and are not involved in permutational 
exchanges, either within the layer itself or with atoms of 
the second layer. Both the first and second layers are 
slightly compressed as further layers form on top. At its 
highest compression, the second layer has a half-width at 



half maximum around 1 A, and corresponds to an effec- 
tive coverage of approximately 0.050 A -2 , i.e., the second 
layer is already liquid. This contradicts the assumption 
of two solid adsorbed layers of helium on the substrate, 
modelled as in this work (with the same choice of parame- 
ters) , made in Ref. d, as the potential ([2]) with the choice 
of parameters made here, only leads to the formation of 
a single solid layer. 



Superfluidity 

Unlike structure and energetics, the superfluid prop- 
erties of a systems studied by computer simulation are 
strongly affected by its finite size, i.e., the number of 
particles. In particular, phase transitions are smeared, 
and one will typically observe a finite value of the su- 
perfluid fraction ps for temperatures above the critical 
temperature T c . 

In this study, given the crudeness of the model uti- 
lized we are not interested in providing accurate values 
of the transition temperature, but rather to assess reli- 
ably whether an adsorbed film of a given coverage can be 
expected to display superfluidity below a reasonably well- 
defined temperature. Because of the two-dimensional 
(2D) nature of the system studied here, the superfluid 
transition is expected to fall in the Kosterlitz-Thouless 
universality class [l8j], with a sharp universal "step" at 
T c . Finite-size effects affect mostly estimates obtained 
in the vicinity of T c , whereas for T < T c , the results for 
a finite system are quantitatively representative of the 
thermodynamic limit, where ps ~ Ps{T — 0) [unity for 
a translationally invariant system] . 

The most recent estimate for lf D is 0.655 K, based 
on the Aziz pair potential utilized here, which affords 
fairly closed agreement with experiment for bulk liquid 
4 He [l(J [lH. Computer simulations of an adsorbed 4 Hc 
monolayer on Lithium 0], yielded evidence that such a 
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FIG. 6: Snapshot of an instantaneous many-atom configura- 
tion for a 4 He film of coverage 0=0.299 A -2 at T=0.5 K. 
Shown are particle world lines projected on the xz plane. 
Clearly identifiable are only the two bottom layers, the sec- 
ond being liquid. Quantum many-particle exchanges involv- 
ing atoms in the second layer are largely restricted to that 
layer. 



film essentially mimics the physics of 2D 4 He, and be- 
cause the substrate considered here is considerably more 
attractive than Lithium, we can expect the same to be 
true here. Indeed, the superfiuid fraction of the equilib- 
rium stable monolayer film at T=0.5 K, is close to 100%, 
consistently with its behaviour being fairly close to that 
of purely two-dimensional 4 He (such a monolayer is con- 
siderably narrower than on a Li substrate, see Ref. 0). 

As mentioned above, as the coverage is increased the 
first layer crystallizes, and superfluidity disappears, to 
reappear only when a stable liquid second layer forms, at 
a coverage 9 2 l- As shown in Fig. El the superfiuid frac- 
tion at 6 = 9 2 l at T=0.5 K is close to 33%, i.e., consistent 
with an underlying non-superfluid crystalline layer, and 
an entirely superfiuid second layer on top. That particles 
in the first layer are not involved in quantum-mechanical 
exchanges is inferred by histograms of permutation cy- 
cles, as well as from the lack of overlap between the first 
and second layer, as shown in Figs. [3]and[6l 
Given that no evidence is observed of crystallization for 
layers above the second, and that the superfiuid transi- 
tion temperature is expected to be generally an increasing 
function of the film thickness (the superfiuid transition 
temperature should approach the bulk value T=2.177 K 
in the thick film limit), it seems reasonable to posit that 
T=0.5 K should be sufficiently low to render finite-size 
effects essentially unnoticeable. 

At coverages above 9 2 l, quantum exchanges involving 
particles in the same layer (other than the first) occur at 
T=l K already, interlayer exchanges being more frequent 
between the third and fourth layer, which overlap most 
substantially (as shown in Fig. [5]). Consequently, the 
superfiuid fraction is finite for any stable film consisting 
of more than one layer, confirming that all layers above 
the first are liquid. Results for the 4 He superfiuid frac- 



tion ps at T—0.5 K are shown in Fig. [5l The behaviour 
as a function of 9, for the coverages above 9 2 l is oscil- 
latory, with local minima indicating layer "completion" . 
However, ps stays finite at all coverages. As the coverage 
increases, the superfiuid fraction is expected to approach 
100% in the thick film limit. 



Single He impurity. 

The physics of a single 3 He atom dissolved in an ad- 
sorbed 4 He film is of interest for a number of reasons. 
First of all, He is never isotopically pure, as contami- 
nating 3 He impurities are always present and their effect 
on the behaviour of 4 He in confinement, or restricted 
geometries, is not completely understood. Moreover, iso- 
topic liquid helium mixtures are of considerable funda- 
mental interest, as one of the simplest Fermi-Bose mix- 
tures accessible in nature. Their phase diagram has been 
extensively investigated experimentally, including in the 
confines of porous glasses such as aerogel 

The most basic question that one can ask, is in what 
region of an adsorbed 4 He film of varying thickness will 
a 3 He atom position itself. The wave function of a sin- 
gle 3 He atom embedded in an adsorbed 4 He monolayer 
is expected to be more "spread out", i.e., the atom is ex- 
pected spend more time further away from the substrate 
than a 4 He, because of its lighter mass. In the case of 
a few adsorbed 4 He layers, it is presumed that the 3 He 
atom will reside in the top layer, near the free surface of 
the film, where the three-dimensional density approaches 
that of bulk liquid 4 He. Conceivably, the 3 He atom might 
even "float" on top of the outmost 4 He layer. 

In order to assess these issues quantitatively, we have 
carried out simulations of a single He impurity at the 
same helium coverages considered for the case of a pure 
4 He film. At all coverages explored here, the 3 He atom 
is found to reside on the top layer of the film. Fig. [7] 
shows the probability density of position in the direction 
perpendicular to the substrate, for a single 3 He atom em- 
bedded in adsorbed 4 Hc films of various coverages, cor- 
responding to one, two and three adsorbed layers. 
In all cases shown, the 3 He and 4 He density profiles over- 
lap substantially, to indicate that the 3 He impurity is 
essentially embedded in the 4 He layer; it does not, e.g., 
float atop. On average, the lighter atom lays only very 
slightly further away from the substrate than a 4 He atom 
on the same layer, although this displacement becomes 
more pronounced as the number of adlayers increases. 
This is attributable to the weakening of the attraction 
exerted by the substrate on the 3 He atom (laying near 
the surface) as the film thickens, and the enhanced quan- 
tum delocalization of 3 He with respect to atoms of the 
heavier isotope. In the limit of a thick 4 He film, the 3 He 
impurity is expected to remain near the surface, where 
the presence of a bound state has been predicted theo- 
retically 2(], 21, 22 1 and verified experimentally 23, 24 1. 
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for different values of coverage. 

It is found that only the first adlayer is solid, whereas 
the successive ones are liquid. The solid "inert" layer 
layer is fairly well separated from the second, atoms are 
highly localized and no quantum exchanges (both intra- 
and inter-layer) take place. On the other hand, successive 
layers are liquid-like, with no sharp separation between 
layers above the third. In the second layers, quantum ex- 
changes are mainly confined to the layer itself, whereas 
interlayer exchanges occur between third and fourth lay- 
ers. 

A single 3 He impurity is always observed on the top layer, 
not near the substrate. This appears to be in contrast 
with the conclusion of Ref. where the existence of 
a homogeneous 3 He- 4 He mixture was predicted near the 
substrate, based on a density functional calculation. On 
the other hand, the results shown here lend credence to 
the scenario proposed in Ref. [25|, according to which 
the heavier isotope in the mixture, 4 He, would be preva- 
lently found in the vicinity of the surface of, e.g., aerogel 
strands. The wave function of the 3 He atom, however, 
is not significantly more spread out than that of a 4 He 
atom on the same layer, i.e., the 3 He atom does not spend 
a significantly fraction of time at greater distances from 
the substrate. 
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FIG. 7: Probability density of position in the direction per- 
pendicular to the substrate, for a 3 He atom embedded in a 
4 He film of coverage 0=0.100 A~ 2 (top, single layer), 0.178 
A~ 2 (middle, two layers) and 0.244 A~ 2 (bottom, three lay- 
ers). Dashed lines show the corresponding probability density 
for a 4 He atom in the same layer where the 3 He atom resides, 
which is always the top one. Units on the ordinate axis are 
arbitrary. 

CONCLUSIONS 

We have studied by accurate, first principles Quantum 
Monte Carlo simulations the physics of an adsorbed film 
of 4 He over a model glass substrate, at low temperature. 
We utilized the most commonly adopted model of a glass 
substrate, neglecting the disorder and irregularity of its 
surface. We considered films up to four layers thick, and 
also looked at a single 3 He impurity diluted in the film, 
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